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Electron energy distribution in a mesoscopic AC-driven diffusive wire generally is not character- 
ized by an effective temperature. At iow temperatures, the distribution has a form of a multi-step 
staircase, with the step width equal to the field energy quantum. Analytic results for the field fre- 
quency high and low compared to Thouless energy are presented, while the intermediate frequency 
regime is analyzed numerically. Manifestations in the tunneling spectroscopy and noise measure- 
ments are discussed. 



In mesoscopic micron-size wires, at low temperatures, 
the electron-phonon energy relaxation time T e - P h can ex- 
ceed the time td of diffusion across the wire [1]. In this 
regime, when electron cooling is primarily diffusive, the 
system temperature has a very short thermal response 
time to microwave radiation, determined mainly by the 
diffusion time T£> . This makes such diffusion-cooled wires 
an excellent system for fast bolometric detection of tera- 
hertz radiation [2,3]. 

When the wire dimension, or electron temperature T e , 
decreases further, one arrives at the situation, demon- 
strated in Saclay experiments [4,5], when the electron- 
electron energy relaxation is longer than the diffusion 
time td. When this happens, the electron dynamics can 
be considered as purely elastic. Importantly, the effect 
of microwave field on electron energy distribution in this 
case is not described by an effective temperature. Since 
the energy of the field is absorbed in discrete quanta Hu>, 
and there is not enough time to redistribute it between 
electrons while they move around in the wire, one expects 
a staircase-like energy distribution to emerge, consisting 
of htu wide steps. The steps will be pronounced when the 
temperature of electrons is below huj. 

In this article, our goal is to develop a general approach 
that allows to analyze such a situation. We shall focus 
on the regime when the field frequency is small compared 
to the elastic scattering rate due to disorder, which cor- 
responds to the experimental situation [2,3,5]. Since in 
this case the transport is diffusive on the field oscilla- 
tion time scale 2n/u), the system can be described by 
momentum- averaged Greens functions, similar to Usadel 
theory of disordered superconductors. We develop a gen- 
eral framework to analyze the energy distribution, using 
Keldysh Greens functions, and then apply it to obtain 
analytic results in the two regimes, when the field fre- 
quency is small and large compared to Thouless energy 
Et = H/td- In the intermediate regime, ujtd ~ 1, we 
present numerical results. 

The energy distribution of this form can be probed by 
tunneling spectroscopy, similar to the DC transport sit- 
uation, when a double step structure arises from mixing 



of the lead Fermi functions at nonequal chemical poten- 
tials [4] . Another manifestation we consider is in the shot 
noise, where the structure of the energy distribution af- 
fects the noise power. 

Before discussing the AC mesoscopic transport prob- 
lem, let us recall, for comparison, the basic facts about 
the DC transport. At low temperatures, when electron 
energy relaxation is slow, charge transport is mainly con- 
trolled by elastic scattering due to disorder. At high con- 
ductance, when the localization effects are negligible, the 
sample properties can be fully described by a scatter- 
ing matrix. All the statistics of transport in this case are 
controlled by a single parameter, the sample conductance 
(see review [6]). 

In the DC regime, the energy gained by an electron 
moving in a stationary electric field depends only on the 
total electron position displacement and is independent 
of the trajectory shape. As a result, the electron energy 
distribution in a wire carrying a DC current does not de- 
pend on microscopic parameters and is fully determined 
by the external voltage and the wire geometry. The en- 
ergy distribution in the DC case was explored by Na- 
gaev [7,8], who obtained a position-dependent mixture of 
two Fermi functions in the limit of slow energy relaxation, 
and a single Fermi distribution with position-dependent 
effective temperature in the limit of fast relaxation, and 
studied the effect on the shot noise. More recently, 
Nagaev [9] has shown that the time-dependent current 
fluctuations of all orders can be expressed through the 
noncquilibrium electron energy distribution. 

In the AC case, in contrast, the energy gained by an 
electron depends on its whole trajectory. The AC trans- 
port can thus reveal additional information about phase- 
sensitive effects in electron dynamics. It was pointed out 
by Lesovik and Levitov [10] that even a slowly varying 
AC field leads to a new transport effect, photon-assisted 
noise, observed by Schoelkopf et al. [11] and recently by 
Glattli et al. [12]. The photon-assisted effects can be ex- 
pected to become more drastic and interesting when the 
field frequency is comparable to Thouless energy. 

In this paper, a semiclassical treatment of ac trans- 
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port in mesoscopic wire is adopted. In a time-dependent 
external field, the uncertainty principle restricts the abil- 
ity to measure the energy of an electron as a function of 
time, and the time-dependent energy distribution func- 
tion should be carefully defined. In this work, we employ 
the Wigner distribution function and show that in the 
presence of an AC field, the time-averaged electron dis- 
tribution consists of hu> steps. Subsequently, we discuss 
the manifestation of such distribution in tunneling spec- 
troscopy [4] and in noise. 

Let us now turn to the analysis of the energy dis- 
tribution. Non-equilibrium electrons in a diffusive con- 
ductor can be fully described by the retarded and ad- 
vanced Greens function G R (t, t', r, r') and G A (t, t', r, r'), 
and Keldysh function F(t, t',r, r') defined in [14]. In a 
diffusive conductor, many quantities can be calculated 
from momentum-averaged Greens function (i.e., taken at 
equal points in space, r = r'). Then, since for nonin- 
teracting electrons, the functions G R (t,t') and G A (t,t') 
are independent of the external fields, any quantity of in- 
terest can be expressed entirely in terms of the Keldysh 
function F(t, t', r). 

The Keldysh function F(t\,t 2 , r) contains information 
about both energy distribution and the single-particle 
density of states. To separate the density of states, we 
introduce the function /(^i, ^2, r) (see, e.g., [17]), 

F = G R -G A -2G R *f + 2f*G A , (1) 
where / * g stands for the convolution 

oo 

(f*g)(t,t') = J f(t,t")g(t",t')dt'. (2) 

— oc 

Then, the details of the electron spectrum arc hidden 
in the Green's functions G R ' A (t\, t 2 ), while the function 
f(t\,t 2 ,r) depends only on the energy distribution. 

Consider a disordered wire of length L with diffusion 
coefficient D, connected to the leads which serve as a 
source of equilibrium electrons. The wire is subject to 
electric field which we describe as an external vector po- 
tential 



A x (t) = (cU/Lu) cos cut , 



(3) 



where U is the amplitude of voltage across the wire, and 
uj is the frequency of the external field. At frequencies 
low compared to the plasma frequency and Maxwell re- 
laxation rate, the longitudinal component of electric field 
is screened. Thus, the electric field can be described 
through the vector potential A(r, t) only. 

The equal point Keldysh function f{t\,t 2 ,r) satisfies 
the two-time diffusion equation [15-17] 



' d 
dt + 



D 



-iV - -A(h,t 2 ,r) 



/(ti,t 2 ,r) = 0, (4) 



with£ + = (t 1 +t 2 )/2 &ndA(t 1 ,t 2 ,r) = A(*i, r)-A(t 2 , r). 
This equation has to be supplemented with the boundary 
condition at the leads: 



f(t u t 2 ,x = 0) = f(t 1 ,t 2 ,x = L) = f F (h -h 



(5) 



where f F {ti — t 2 ) is the Keldysh function of an equilib- 
rium Fermi gas. 

It is convenient to rewrite the vector potential differ- 
ence as 



A x (h,t 2 ,x) = A X (t + ^j - A x (t-^) 

= h - 1 2 . 



(6) 
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Passing to the Wigner representation, we replace r by 
id € . Now, using the general operator relation 



cxp (a-^j *( e ) = *( e + a ) - 

valid for an arbitrary function <&(e), one arrives at 

icU 



A x (h,t 2 ,x) 



sin(wt)I? [i 



where is a finite difference operator 

^/(e)= /(£ + " /2) - /(£ -" /2) 

UJ 

Next, we perform variable rescaling, 

t — ► tw, x^x/L, e^>e/eU 



(7) 



(8) 



(9) 
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i.e., measure time in the units of tu^ 1 , x in the units of L, 
and energy in the units of eU. After that, the equation 
for the distribution function becomes 



d 



1 



dt ujtd \dx 



d 



— sin W u 



f(t,e,x) = 0, (11) 



where td = L 2 /D is the time of diffusion across the wire. 

For a wire connected to the leads, the boundary con- 
dition to the Eq. (11) takes the form 



f(t, e, x = 0) = f(t, e,x=l) = n F (e) . 



(12) 



where np(e) is a Fermi distribution. For a wire that has 
been in equilibrium with the leads before the field was 
turned on, the initial condition is 



f(t = 0,e,x) = n F (e) 



(13) 



for any x. 

Note that the operator P w in Eq. (11) relates the val- 
ues of the distribution function /(e) only for the energies 
that differ by ±^Huj. Thus, at zero temperature the sin- 
gularity at e = in the boundary condition (12) propa- 
gates to the energies e„ = ^nhu. Also, since the Fermi 
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distribution rip(e) is flat for e < and e > 0, the distri- 
bution function is also flat for e„ < e < e n +i. Thus, the 
profile of f(t, e, x) is a series of steps at energies nhuj/2. 
In fact, only the even steps (e2n = nfwj) corresponding 
to absorption of individual field quanta, survive after av- 
eraging over time. 

For slow field, luto <C 1, one may neglect the time 
derivative in Eq. (11). The solution to Eq. (11) is then 

f(t,e,x) = [(l-aOexp(a£> u siiit) (14) 
+ xexp {{x — l)T> u sini)] np(e) . 

The exponent of the operator T> w can be found using 
Fourier transform in the energy domain. The result is 

oo 

exp(aX> w )$(e) = ^ $(e - nw/2) J n (2a/w) , (15) 

n— — oo 

where $(e) is an arbitrary function, J n (x) is the Besscl 
function of nth order, and the sum runs over all inte- 
ger n. As we discuss below, one is mostly interested in 
the time-averaged distribution function, because it can 
be directly measured by tunneling spectroscopy. Sub- 
stituting Eq. (15) into Eq. (14), and using the formula 
J2fe(2asint) = J%(a) to average over t, one finds, with 
the original units restored, 

fo{e,x) = (l-j)F (e,x) + jF (e,L-x) , (16) 

where 



F (e,x)= J2 J fe(^)' 



k>e/uj 



Lhuj 



(17) 



(We mention that an expression of the form similar to 
(14) for nonaveraged function / was used by Altshulcr et 
al. in the calculation of time-dependent noise [18].) 

To interpret the result (16), note that the work per- 
formed by a slow field on electrons that travel a distance 
x along the wire, is eUx/L. Then, according to [19], the 
probability of absorbing n field quanta is J 2 (u)- The fac- 
tor (1 — x/L) in Eq. (16) is the number of the electrons 
at cross-section x coming from the left lead (cf. [7]). The 
second term in Eq. (16) describes electrons that come 
from the right lead. 

In the limit u> <C ell , one may use the asymptotic form 
of the Bessel function at x/u <~ n » 1, which gives 



F (e,x) = 



1, 
l 

rr 

o, 



- cos 1 e, 



e < -1 
|e| < 1 
e> 1 



e = 



Le 

xeU ' 



(18) 



This result can also be derived by time averaging the two- 
step distribution found in [7] over the time-dependent 
voltage difference. 

For fast field, uitd S> 1, the first term in Eq. (11) is 
the most important. Thus, the electron distribution is 



almost time- independent. (The time-dependent part of 
f(t,e,x) is proportional to (wtd) .) Projecting out the 
time-dependent part of the distribution function, i.e., av- 
eraging the Eq. (11) over field period, one finds the equa- 
tion for the time-averaged electron distribution: 



d 2 

L _X> 

dx 1 2 » 



f(e,x)=0. 



(19) 



Using Fourier transform with respect to energy, one can 
find the solution of Eq. (19) in a closed integral form. It 
is more instructive, however, to solve this equation in the 
limit hu> <C ell, where one can replace by d e . This 
brings Eq. (19) to the form of a Laplace's equation in the 
two-dimensional strip < x < 1, — oo<e<oo. Using 
the function 



w = cxp(7ri(a; + iV2e)) . 



(20) 



one can conformally map this strip onto half-plane 
Imw > 0. The boundary condition (12) on the line 
Imw = is, for T e = 0: 



foo{w) 



0, \Rcw\ < 1 , 

1, | Rew\ > 1 ' 



(21) 



This boundary value problem is solved by an imaginary 
part of the analytic function 



/ooH = Im- log ( \— - 

7T V 1 + W 



(22) 



Restoring the original dimensional units, one obtains 

/ oo( e,, ) = lcotan-^ Sinh ^;/f) '| (23) 

Unlike the low frequency solution (16) with .Fb^i x) given 
by Eq. (18), the high frequency distribution Eq. (23) is 
non-zero at large e. Physically, this difference is due to 
tha fact that maximal energy that can be gained from the 
field is given by the bias amplitude eU only in the DC 
regime. In contrast, in a varying field, there is a finite 
probability that electron diffuses several times back and 
forth in phase with the field, thus gaining energy that 
exceeds eU. 

Eqs. (16), (18) and (23) are correct in the limit ui <C 
eU. For finite values of uj/eU these results approximate 
the average profile of the step-like electron distribution. 
To obtain the distribution at intermediate frequencies, 
ujtd — 1, Eq. (11) was solved numerically for different 
values of eU/uo. The results are shown on Fig. 1. From 
that figure one may see that a crossover from the low- 
frequency behavior to the high-frequency behavior occurs 
at lotd <~ 100. To understand the origin of this large 
numerical factor, note that the electron distribution re- 
laxes at t — > oo, according to the diffusion equation, as 
exp(— lit), where /x is the lowest non-zero eigenvalue of 
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a diffusion operator: /i = tt 2 /to- Assuming that the 
crossover occurs when the relaxation time is of order of 
field period, 2tt/lu, one finds lute, ~ 2-7T 3 w 60, in a qual- 
itative agreement with the numerical data. 
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FIG. 1. The energy distribution at the wire midpoint 
x = L/2, calculated numerically for huj = 0.2eU, for different 
values of lutd- The height of the Tiuj steps is a function of 
external voltage and other system parameters. The envelope 
of the staircase structure changes as the frequency of the field 
increases, demonstrating a crossover from the slow to the fast 
field regime. 

To characterize the electron distribution, one may com- 
pute its second moment, or shot noise. For low frequen- 
cies, the answer can be derived just by averaging over 
time the well known result for the shot noise [7]: 



S = -eG\U(t)\ = —eGU 



0.424eG(7 . 



(24) 



where G is a conductance of the wire. For high frequen- 
cies, one may neglect the time dependence of the distri- 
bution function and compute the noise as [7] 



AeGU 



de /(»(£, x)(l - /oo(e,a;)). (25) 



-oo 



Evaluating the integral, one arrives at 
2lV2 



2tt 3 



■C(3)eGt7 « 0.575 eGU , 



(26) 



where £(3) is Riemann zeta-function. The value given 
by Eq. (26) is higher than that in Eq. (24), because of the 
contributions of high energy tails of the distribution (23). 

The electron distribution can also be probed by super- 
imposing a DC voltage U^ c , which splits every step into 
two, separated by elide- When this splitting becomes 
equal to hu), the "resonance" between two steps leads to a 
singularity in the shot noise and higher current moments, 
predicted for short wires in [10] and observed in [11]. 
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FIG. 2. Schematic representation of the measurement [4], 
in which the energy distribution in a mesoscopic wire (W) is 
probed by a NS tunneling junction between the wire and the 
superconducting probe (SC). By changing voltage V across 
the junction, one may scan the energy of tunneling electrons. 
AC voltage U(t) applied between the leads [L and R), creates 
a nonequilibrium energy distribution in the wire. 

Finally, we discuss how the structure in the energy 
distribution considered above will manifest itself in the 
tunneling spectroscopy measurement [4,5]. In this setup, 
a short segment of the sample at point i"o is connected 
to a metallic probe by a tunneling junction (see Fig. 2). 
The tunneling density of states of the probe must con- 
tain a sharp feature that serves as a "pointer" , allowing 
to probe electron distribution. (In [4], the probe was 
a superconductor, and the BCS square root singularity 
near the superconducting gap was used.) Applying the 
probing voltage V to the junction, allows one to scan the 
energy domain. The electron energy distribution in the 
sample is then extracted from the I — V curve of the 
tunnel junction. 

To derive the relation between the tunneling current 
and the energy distribution in a nonequilibrium system, 
one may use the standard tunneling Hamiltonian 



P ,k 



(t)6k(t) +h.c. 



(27) 
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(see, e.g., [13]). Here a p (t) and b^(t) are electron opera- 
tors of the sample and the probe. The tunneling current 
operator in this formalism is 



ir(t) = i^T Pjk a+(i)6 k (i) + h.c. . 



(28) 



The average tunneling current I{t) can be found from 
Kubo formula: 



lit) 



dt' 



H T (t'),i T (t) 



(29) 



Using the Kcldysh approach [14], one may express the 
tunneling current I(t) in terms of retarded and ad- 
vanced Green's function G R ( A ^ (t\,t 2 ) and Keldysh func- 
tion F{t\,t 2 ) defined in [14]: 



J(t) oc Re [F s * G R — F p * G 



s J lt'=t 



(30) 



(subscripts s and p denote the sample and the probe, all 
functions are taken at space point ro, where the tunnel- 
ing occurs). The star denotes the convolution defined by 
Eq. (2). 

It is now convenient to express all quantities in fre- 
quency domain. In the nonintcracting system, the func- 
tions G R ' A {t\, t 2 ) depend only on the difference of their 
arguments ti — t 2 . Their Fourier components are related 
through the Kramcrs-Kronig relation to the density of 
states v(uj) defined by 



TTv{uu)=lmG R {uj) = -luvG A {uo). 



(31) 



The function f(ti,t 2 ) in general depends on both argu- 
ments. Wc define its Fourier transform as 



(32) 



Thus, fu{e) is a Fourier transform of the Wigner distri- 
bution 



/M) = f fit 



+ T/2,t-r/2)e l€T dr. 



(33) 



Then, the Fourier components of the tunneling cur- 
rent (30) are given by 



lui oc 



oc 

J deu s (e + |) v v (e + |) [/ a , w ( e ) - / p , w (e)] . (34) 



In an equilibrium system, the function / w (e) coincides 
with the Fermi distribution [14]: f u {e) = 2nn F (e)5(uj). 
Thus, /(e, t) can be identified as the conventional energy 
distribution function. 

The probing voltage V simply shifts energy levels of 
the probe by eV. Thus, the tunneling current is: 



T w {eV) oc J de v s (e + ^ u p (e + ^ - 
x [/ a , w (e)- / PlW (e + eV)] . 



eV) 



(35) 



Eq. (35) shows that in general, the relation between 
the tunneling current and energy distribution function is 
not straightforward, because all factors in (35) depend on 
the current frequency u. (One cannot extract the energy 
distribution from tunneling measurements if both v s (€) 
and Vp{e) are flat.) The w-dependence of the product 
v s (e + uj /2)is p (e + u /2 + eV) leads to a non-local relation 
between energy distribution and tunneling current. This 
non-locality does not allow to measure Wigner function 
directly. However, for u> = 0, the average tunneling cur- 
rent is proportional to the average electron distribution: 

I(eV) ex J deu s {e)v p {e + eV) [f s (e) - f p (e + eV)} , 

(36) 

where the bar denotes averaging over one field period. 
Thus, the tunneling spectroscopy technique of [4] allows 
to measure directly the time-averaged electron energy 
distribution. 

The effect of finite temperature can be understood as 
follows. At finite electron temperature T e , the singularity 
in the boundary condition (12) is smeared. The result- 
ing smearing of the step-like structure is negligible for 
fiuj > T e . For a wire length L = l/im, the electron mean 
free path I = lOnm, and Fermi velocity vf = 10 7 m/s 
the diffusion time is td = L 2 /vpl = 10 -11 s. The fast 
field regime occurs for luto > 100, i.e., for the frequen- 
cies / = uj/2tt > 1 THz. The optimal voltage amplitude 
is of the order of IOHuj, i.e. V ~ 50 mV. The temper- 
ature must be less than fkv/kB, i.e., T < 10 K. These 
conditions appear realistic, and can be fulfilled in an ex- 
periment. 

The simple theory presented in this paper does not take 
into account the effect of non-equilibrium density fluctua- 
tions in the wire on the density of states of the supercon- 
ducting probe through long-range Coulomb forces [20]. 
The treatment of this effect for an AC field is beyond the 
scope of this paper. However, since the non-equilibrium 
fluctuations are also due to periodic field, the resulting 
tunneling current will also have singularities at quantiza- 
tion energies nhui, although the shape of these singular- 
ities may differ from simple steps. 

To summarize, the energy distribution of a mesoscopic 
AC-driven wire is not described by an effective tempera- 
ture. The most interesting feature of the electron distri- 
bution is the hco steps due to quantization of absorbed 
energy. Those steps are present at the external field 
frequencies larger than electron temperature. The cal- 
culated envelope in the fast-field regime distribution is 
qualitatively different from that of the slow-field regime. 
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Experimental manifestations of these effects experimen- 
tally are discussed. 
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